Pair counting, pion-exchange forces, and the structure of light 

nuclei 

R. B. WiringaQ 

Physics Division, Argonne National Laboratory, Argonne, Illinois 60439 

(Dated: February 9, 2008) 

Abstract 

A simple but useful guide for understanding the structure of light nuclei is presented. It is 
based on counting the number of interacting pairs in different spin-isospin (5, T) states for a given 
spatial symmetry, and estimating the overall binding according to the sum of ai-ajTi-Tj expectation 
values, as suggested by one-pion-exchange. Applied to s- and p-shell nuclei, this simple picture 
accounts for the relative stability of nuclei as A increases and as T changes across isobars, the 
saturation of nuclear binding in the p-shell, and the tendency to form d, t, or a subclusters there. 
With allowance for pairwise tensor and spin-orbit forces, which are also generated or boosted by 
pion-exchange, the model explains why mixing of different spatial symmetries in ground states 
increases as T increases across isobars, and why for states of the same spatial symmetry, the ones 
with greater S are lower in the spectrum. The ordering of some sd-shell intruder levels can also 
be understood. The success of this simple model supports the idea that one-pion-exchange is the 
dominant force controlling the structure of light nuclei. 

PACS numbers: PACS numbers: 21.10.-k, 21.45. +v, 21.60.-n 
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I. INTRODUCTION 



The last decade has seen significant progress both in the characterization of realistic two- 
and three-nucleon interactions, and in the ability to make accurate many-body calculations 
with these models. Nucleon-nucleon potentials such as Argonne vx$ |l| , CD-Bonn , and 
the Nijmegen models Q| reproduce AW scattering data extremely well; when combined with 
three-nucleon forces such as the Illinois or Tucson- Melbourne potentials 0, |(| and accurate 
many-body techniques, nuclear binding energies up to A = 12 can be reproduced. Seven 



different many-body methods are in superb agreement for the binding ener gy o 



realistic AW force jzj], while Green's function Monte Carlo (GFMC) 



shell model (NCSM) [12, ll3L Il4j|. and coupled-cluster methods (CCM) |l5l llq are making 




4 He with a 



111 ], no-core 



very successful ab initio calculations for p-shell nuclei. This progress allows us to study the 
interplay between nuclear forces and nuclear structure in an unprecedented manner. 

In a recent letter Q| we constructed a series of increasingly realistic force models and used 
GFMC calculations to evaluate the consequences for nuclear structure. This study showed 
that a simple central potential, with the canonical intermediate-range attraction and short- 
range repulsion indicated by S"-wave AW phase shifts, could approximately reproduce the 
triton and alpha binding energies, but failed to saturate in the p-shell, producing stable 
5 He and greatly overbinding the A=6,7,8 nuclei. To obtain unstable 5 He and the general 
saturation of nuclear forces that is evident in the p-shell, it is necessary to have a state- 
dependent force, i.e., one that is attractive in L=even partial waves and repulsive in L=odd 
partial waves. Indeed, state-dependence appears to be more important for nuclear saturation 
than either the repulsive core or the finite range of nuclear forces. To obtain the further 
refinement that 8 Be is unstable against breakup into two alphas requires the addition of a 
tensor force, while the stability of 6 ' 7 Li suggests spin-orbit terms are also needed. Together, 
these are the major operator components required in a realistic interaction to fit S- and 
P-wave AW data. 

Pion-exchange forces play a very significant role here. The spin-isospin dependence of one- 
pion exchange (OPE) makes it attractive in L=even partial waves and repulsive in L=odd 
partial waves, just as required to bind the s-shell nuclei but saturate quickly in the p-shell. 
OPE also is the major source of the tensor force, and iterated tensor interactions between 
three or more nucleons provide a large enhancement to spin-orbit splitting in nuclei jlj. 
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In GFMC calculations of A < 12 nuclei with realistic interactions, the expectation value of 
the OPE potential is typically 70-75% of the total potential energy ^. The importance 
of pion-exchange forces is even greater when one considers that much of the intermediate- 
range attraction in the NN interaction can be attributed to uncorrelated two-pion exchange 
with the excitation of intermediate A(1232) resonances [3]. In addition, two-pion exchange 
between three nucleons is the leading term in 3N interactions, which are required to get the 
empirical binding in light nuclei j^. In particular, the 3N forces provide the extra binding 
required to stabilize the Borromean nuclei 6 ' 8 He and 9 Be. 

The thesis of this paper is that by counting the number of different spin-isospin (S, T) 
pairs that occur in a given nuclear state of specific spatial symmetry, and multiplying by a 
numeric strength taken from the OPE operator ovOj'Tj-Tj, a very good measure of the binding 
energy is obtained. This works both for the relative energy between different states in the 
same nucleus, and between different nuclei. The idea is akin to the supermultiplet theory 
of Wigner |20j which focused on the symmetry aspects of light nuclei, but assumed forces 
that were primarily central and space-exchange in character. The present study benefits 
from the extensive recent progress in fully realistic calculations mentioned above. There is 
also common ground with the recent work by Otsuka and collaborators (3, 2 within the 
framework of traditional shell model that emphasizes the importance of OPE spin-isospin 
and tensor interactions in determining how single-particle energy levels shift as shells are 
filled. 

This simple guide, supplemented by our knowledge of tensor, spin-orbit, and Coulomb 
forces, describes the general structure of light nuclei in considerable detail. The model 
explains the growing binding as A increases, the saturation of binding going from the s- 
shell to the p-shell, the relative stability as T varies across isobars, and the tendency to 
form d, t, and a subclusters in the light nuclei. It explains why mixing of different spatial 
symmetries in ground states increases as T increases and why for states of the same spatial 
symmetry, the ones of higher 5* are lower in the spectrum. The same logic can also be used 
to understand the ordering of some sd-shell intruder levels in some of these nuclei. 
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II. PAIR COUNTING 



The total number of pairs in a nucleus, P A = A(A — l)/2, can be subdivided into pairs of 
specific spin and isospin P A (ST) where S = or 1 and T = or 1. Starting with the square 
of the expression J2i T i/2 — T A , where T A is the total isospin of the nucleus, and using the 
projection operators (1 — Ti ■ Tj)/A and (3 + 73 • Tj)/A for T = and 1 pairs, respectively, the 
total number of such pairs of given isospin in a nucleus can be shown to depend only on A 



and Ta 



23j: 



P A (10) + Pa(00) = h A 2 + 2A- AT A (T A + 1) ] , (1) 

o 

P A (ll) + Pa(01) = h 3A 2 - 6A + 4T A (T A + 1) ] . (2) 

o 

A similar pair of equations can be obtained for the total number of S = or 1 pairs in terms 
of the total nuclear spin Sa, assuming spin is conserved, i.e., before configuration mixing by 
tensor forces and correlations. An additional expression can be obtained for the difference, 
P[ n ], between the number of symmetric (even) and antisymmetric (odd) pairs for a given 
spatial symmetry state specified by the Young diagram [n]: 

P A (10) + P A (Q1) - P A (ll) - Pa(00) = P[n] . (3) 

For example a [3] symmetry state has three symmetric pairs, p™ = 3; a [111] state has three 
antisymmetric pairs, P[m] = —3; and a [21] state has one symmetric, one antisymmetric, 
and one mixed symmetry pair (which does not contribute here), giving pt 2 i] = 0. Together 
we have four independent relations for four unknowns, which can be rearranged to give: 

P A (11) = ^[2P A - P[n] -^A + S A (S A + 1)+T A (T A + 1)] , (4) 

P A (10) = \[ Pa+P[u] +S a (S a + 1)-T a (T a + 1)} , (5) 

Pa(01) = \[ P A +P [n] - S A (S A + 1) + T A (T A + 1) ] , (6) 

Pa(00) = - A [ - P[n] + ^A-S A (S A + 1)-T A (T A + 1)} . (7) 

The simple energy measure we propose using is obtained by multiplying the number of 
pairs of each type with the expectation value of the spin-isospin operator • ajT- t ■ Tj coming 
from one pion-exchange: 

Eqpe = C [ P A (11) - 3P A (10) - 3P A (01) + 9P A (00) ] , (8) 



where C is a constant in MeV. A value C ~ 1.5 MeV gives a reasonably good average scale 
factor. This expression reflects the fact that S'-wave AW interactions are attractive while 
P-wave interactions are repulsive. It does not attempt to differentiate between 1 Sq and 3 Si- 
3 Di channels, when in reality the former is just unbound and the latter produces a bound 
deuteron, thanks largely to the OPE tensor force. However it does reflect the large difference 
between the weakly repulsive 3 Pj channels and the strongly repulsive 1 P\ interaction. We 
will find that this simple expression does a remarkably good job of predicting overall trends 
in binding and relative stability for s- and p-shell nuclei, as well as explaining a variety of 
observed features in the excitation spectra. 

III. ENERGY SPECTRA 

The Pa{ST) and Eope for A=2-5 nuclei are shown in Table HJ (In the following tables, 
we will only show the most neutron-rich member of any isobaric multiplet, e.g. 3 H but not 
3 He, and 6 He but not 6 Be or the isobaric analog states in 6 Li; they should be understood to 
be essentially the same for nuclear forces and differ primarily by the Coulomb energy.) The 
deuteron, 2 H, has only one ST =10 pair, which we assign the strength — 3C. In our simple 
model, the ST =01 dineutron would also be bound, whereas in reality it is just unbound, 
and is not shown. The triton, 3 H, has three pairs, equally divided between ST =10 and 01 
according to Eqs.(jUEj), and thus gets a strength of — 9C, while the alpha, 4 He, has six such 
pairs with total strength of — 18C. 

If we use Eqpe to judge the relative binding of these nuclei, then the d : t : a energies 
should be in the ratio 1:3:6, whereas in reality they are more like 1 : 4 : 13. Of course the 
binding is the result of a cancellation between kinetic and potential energies, and Eqpe is 
essentially a potential measure. In fact, GFMC calculations for the AV18/IL2 Hamiltonian 
give the expectation values for the two-body potential to be in the ratio 1 : 2.7 : 6.3 for these 
nuclei 0], reasonably close to E OPE . The E OPE will be a useful gauge for binding energies 
only if there is something like a virial theorem for nuclei that says the kinetic and potential 
energies are proportional to each other. Fortunately there does seem to be such a relation, 
at least in the light p-shell nuclei, as shown by the results of GFMC calculations displayed in 
Table ITT1 For 6 < A < 12 nuclei, the ratio, Rkv, of kinetic to potential energy expectation 
values varies only from 0.76 to 0.81, while there is a much greater range of 0.75 to 0.90 for 
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the s-shell nuclei. We note that the lowest ratios occur for the most spatially-symmetric 
nuclei 4 He, 8 Be and 12 C as one might expect. 

A further complication is that when realistic tensor forces are included, some fraction of 
the ST =01 pairs will be converted to 11 pairs, and a small fraction of ST =10 pairs to 00 
pairs, due to multi-body correlations. Variational Monte Carlo (VMC) calculations for 4 He 

2J. For 



found the actual distribution of 11 : 10 : 01 : 00 pairs to be 0.47 : 2.53 : 2.99 : 0.01 
our simple model we will focus on the distribution of pairs before such mixing takes place. 

While this very qualitative model is not particularly useful for the s-shell, it starts to 
have some utility in the p-shell. The ten pairs of nucleons in 5 He are divided into six pairs 
within the s-shell, designated "ss" in Table EJ and four pairs with one nucleon in the s-shell 
and one in the p-shell, designated "sp". The ss pairs are distributed exactly as in 4 He, 
while the sp pairs come in just the right combination to give no additional contribution 
to Eqpe- Thus the prediction is that 5 He should have the same binding as 4 He, when in 
fact it is unstable against breakup by ~1 MeV. For cases like this, our simple measure is 
not sufficient to determine stability, but only to indicate a situation that could go either 
way, depending on, for example, how much the virial ratio, Rrv, of Table ITT1 varies. In the 
case of 5 He, the residual attraction from shorter-range NN and 3N forces is not enough to 
overcome the additional kinetic energy that is generated by the requirement of putting the 
fifth nucleon in a p-shell orbital. 

For A=6 nuclei, shown in Table UTT] the s-shell core remains the same, while the number of 
sp pairs doubles, but still with no net contribution to Eope- Effectively, the s- and p-shells 
decouple from each other at the OPE level, which may help us understand why the use of 
an inert core in the traditional shell model approach is valid. The final energy of 6 Li and 
6 He then depends on the last pair of nucleons which is wholly within the p-shell, designated 
"pp" in Table ITTTl Starting with A=6, there are multiple ways of adding up orbital and spin 
angular momenta to get the total J n ; T of a given nuclear state j^; we label them by their 
LS coupling and spatial symmetry, 2S+l L[n], and list all allowed L values. For 6 Li this last 
pair can be part of a 3 S'[42] or 3 D[42] state (essentially a deuteron with orbital momentum 
of or 2 around an alpha core) with an associated strength of — 3C, or the last pair can be 
part of a X P[411] state, which contributes +9C to our binding measure. The total Eqpe for 
6 Li is the sum of the ss and pp pairs, or — 21C for the [42] states and — 9C for the [411] state. 
For 6 He the last pair can be part of either a 1 5'[42] or 1 D[42] state (essentially a spin zero 
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dineutron with L = or 2 around an alpha core) with strength — 3C, or part of a 3 P[411] 
state with strength +C; the corresponding total Eope is —21C or —17C, respectively. 

Thus the prediction of our simple model is that 6 Li and 6 He ground states should have 
[42] symmetry, with about the same energy, and be weakly bound compared to 4 He, which is 
pretty much correct, given the above caveat about unbound dineutron and bound deuteron. 
The experimental spectrum 2^| is shown in Fig. ^ where the levels are labeled by their 
dominant symmetry. The 6 Li ground state is 1.47 MeV below the alpha-deuteron threshold, 
while 6 He is 0.97 MeV below the alpha-dineutron threshold. GFMC calculations indicate 
that much of the binding between clusters is provided by the 3N force; if only the AV18 
NN force is used, 6 Li is stable by 0.6 MeV and 6 He is unstable by 0.3 MeV [lpj. 

Not surprisingly, in the excitation spectrum the D states are higher than the S states, 
because the angular momentum barrier screens some of the overall potential attraction. 



In shell mode 
interaction 1 26 



studies this feature is taken into account by including an L 2 term in the 
, but the spread between different L-states of the same spatial symmetry is 
generally smaller than the spacing between different spatial symmetry groups, and here we 
are after only the most general nuclear structure aspects. Further, the 3 D combination in 
6 Li is split into J=l,2,3 states ordered with maximum J lowest, as dictated by the spin- 
orbit force. The antisymmetric [411] states are several MeV higher in the spectrum and no 
corresponding experimental states have been identified. However our simple model predicts 
that the antisymmetric 1 P[411] state in 6 Li is much higher than the 3 P[411] state in 6 He, 
which suggests that when configuration mixing with tensor forces is done, the admixture of 
these components in the respective ground states will be less for 6 Li than for 6 He. This is 
borne out in VMC diagonalizations with realistic forces where the amplitudes of the different 
components in the ground state are 0.98 : 0.14 : 0.10 for the 3 5'[42] : 3 P[42] : 1 P[411] pieces 
in 6 Li and 0.97 : 0.23 for the 1 S[42] : 3 P[411] pieces in 6 He Q- 

The Pa{ST) and Eope for A = 7 nuclei are given in Table HVl Again we see that the 
s-shell core gives the same contribution as before, and though there are now twelve sp pairs, 
they continue to give no net contribution to Eope- All the action is now in the three pp 
pairs. In 7 Li they can form part of a maximally symmetric 2 P[43] or 2 P[43] state with 
energy contribution — 9C for a total Eqpe = — 27C, which equals the sum of alpha and 
triton energies; experimentally the ground state is 2.47 MeV below this sum, as seen in the 
experimental spectrum of Fig. El Again, much of the binding between clusters is apparently 
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due to the 3N forces; GFMC calculations with AV18 alone produce a 7 Li ground state only 
0.3 MeV below the alpha-triton threshold jlp| . 

As labeled in the figure, the 7 Li states are ordered according to our simple model, with 
2 P[43] and 2 P[43] states lowest, followed by the 4 P[421] states and the start of the 4 P>[421] 
states; VMC calculations confirm that these are by far the dominant components of the 
first eight states. (The lowest five states in 7 Be follow a similar pattern with a Coulomb 
shift; the higher states may not be as well known experimentally.) The 4 P[421] states have 
a net Eqpe — —21C, the same as 6 Li ground state, and they lie just above the threshold for 
breakup into 6 Li+n. The 2 P[421] and 4 P[421] T = 1/2 states in 7 Li have the same spatial 
symmetry but the former contain an admixture of the very repulsive ST=00 pairs, which 
pushes their energy up significantly, and no corresponding experimental states have been 
identified. In contrast, the 2 P[421] T = 3/2 ground state in 7 He (and its isobaric analogs) 
does not have any ST=00 pp pairs, and by our simple model has the same energy as the 
6 He ground state, which is about right. 

The Pa{ST) and Eope for ^4 = 8 nuclei are given in Tables IVl and IVTl In the former we 
show the ss and sp pairs again to remind us that each of the p-shell nuclei has an s-shell core 
contributing — 18C to Eqpe and no contribution from the sp pairs. In 8 Be the six pp pairs 
in the maximally symmetric [44] ground state effectively form a second alpha, so the total 
E OPE is 2 x (-18C) = -36C, compared to -28C for 8 Li and -24C for 8 He ground states. 
This is a fair representation of the spread in the experimental spectrum, shown in Fig. El 
[22]]. The 8 Be ground state is practically degenerate with the energy of two alphas, 8 Li is 
significantly less bound, but is a little more bound than 7 Li (— 27C), and 8 He is somewhat 
less bound, but below 6 He and 7 He (both — 21C). The increased binding for 8 He is essentially 
due to the completion of a second dineutron pair in its [422] symmetry ground state, which 
is worth an additional — 3C in Eqpe- 

The gap between the [44] and [431] symmetry states in 8 Be has the large value of 12C, 
suggesting little mixing between them, and VMC calculations indicate the first + , 2 + , and 
4 + states are ~ 99% pure symmetry [44]. By comparison, the small energy gap of AC 
between the [422] and [4211] symmetry states in 8 He leads to mixed amplitudes of 0.8 : 0.6 
in its ground state [lOj. In 8 Li there are both triplet and singlet states of symmetry [431], 
but the higher spin states fall lower in the spectrum because they avoid the repulsive 5T=00 
pairs; the same is true for the three different spin states of [4211] symmetry, and in 8 Be for 
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the two different spin states of [422] symmetry. 

The Pa{ST) and Eope in A = 9 nuclei are given in Tables IVTTI and IVTTT1 Adding in the 
contribution from the s-shell core, the Eqpe are — 36C, — 30C, and — 24C for the ground 
states of 9 Be, 9 Li, and 9 He, respectively, which is again a very good approximation to the 
experimental spectrum shown in Fig. 0] [27]. The 9 Be ground state is predicted to have 
the same energy as 8 Be or two alphas: the addition of one nucleon to the ground state 
of 8 Be generates four new pp pairs in 9 Be, but with just the right combination to add no 
additional binding to E OPE . Experimentally 9 Be is bound with respect to the threshold 
for alpha-alpha-neutron breakup by 1.57 MeV, which in turn is 0.10 MeV below the 8 Be+n 
threshold. GFMC calculations indicate that the stability of the last neutron is again due to 
3iV forces: whereas the AV18/IL2 Hamiltonian gets 1.9 ± 0.5 MeV for the binding relative 
to 8 Be, AV18 alone is stable by only 0.1 ± 0.4 MeV Q. 

The 9 Li ground state is predicted to be somewhat more bound with respect to 8 Li (— 30C 
compared to — 28C) and experimentally it is stable by 4.06 MeV. On the other hand, 9 He is 
predicted to be the same energy as 8 He (both — 24C) because the last neutron is unpaired; 
experimentally the lowest natural-parity l/2~ state is unbound by ~ 1.2 MeV. However, 
recent experiments indicate the lowest state in 9 He is an unnatural positive-parity l/2 + 
state just above threshold, and there are also many low-lying positive-parity states in 9 Be, 
starting with a l/2 + state just above the 2a + n threshold. These unnatural-parity states 
can be constructed by putting the last nucleon in an sd-shell orbital outside the p-shell 
core. As with the p-shell, the Eope weight factor is such that there is no net interaction 
between an sd-shell nucleon and the core. In the case of 9 He and 9 Be, our simple model 
suggests that the the long-range part of the NN potential does not care what orbital the 
last nucleon goes into; whether a p-shell or sd-shell orbital is more stable depends on the 
residual shorter-range NN interaction, the 3N interaction, and the kinetic energy cost. 

There is a moderate gap of size 6C between the first and second symmetry states in 9 Be, 
which is not as large as the gap in 8 Be; consequently the low-lying states are mostly [441] 
symmetry with relatively small admixtures of [432] components, but not as pure as the 8 Be 
[44] states jjj. The smaller gap between symmetry states in 9 Li leads to more mixing of the 
[432] and [4311] components there. VMC diagonalizations also continue to show that for 
states of the same spatial symmetry, those with higher spin lie lower in the spectrum, again 
due to a smaller presence of repulsive ST=00 pairs. 
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The Pa(ST) and Eope in A = 10 nuclei are shown in Tables ITXl IXl and IXI1 The Eope 
are -27C, -31C, -39C, and -39C for 10 He, 10 Li, 10 Be, and 10 B, respectively. Thus 10 B 
and 10 Be should have the same binding, which experimentally they do at 64.75 and 64.98 
MeV as seen in Figs. and El They are predicted to be about 4-5 MeV more bound than 
8 Be and 9 Be, but experimentally it is more like 7-8 MeV. On the other hand, the prediction 
for 10 Li is that it should be a little more bound than 9 Li, whereas it is unbound by about 
0.25 MeV. Further, 10 He should be bound by several MeV compared to 8 He and 9 He by 
the completion of another dineutron pair, whereas it is unbound by 1 MeV compared to 
8 He. This could be an indication that jj-coupling is starting to be more appropriate as the 
neutron p-shell is completed 2sj], with this last pair of neutrons being a pi/ 2 pair that joins 
at a noticeably higher energy than the first two dineutrons. 10 Li and 10 He are the only two 
nuclei out of 27 in the p-shell (not counting isobaric analogs) that are falsely predicted to 
be stable by our simple model. 

Comparing the different symmetry states in Table IXTl for 10 B, we see that the [442] com- 
ponents are 12C below the [4411], [433], etc., components, so there is very little admixture 
of the latter into the lowest-lying states. However, in 10 Be the gap between the ground state 
[442] symmetry and the [4411] and [433] components is only 4C, so there is a moderate 
admixture into the ground state [jj. The clearest signature for a state of these next spatial 
symmetries in the p-shell would be a 1 + state in 10 Be, expected at ~ 6 MeV excitation, 
but no such state has been observed. However, unnatural-parity states which involve an sd- 
shell intruder, are now low enough in the spectrum for some of them to be particle-stable; 
discussion of these is deferred to the next section. 

In full GFMC calculations, the 3N force starts to make an especially large impact by 
A = 10 in that it starts to reorder some of the states from our simple expectations. Naively 
we would expect the 10 B spectrum to be something like 6 Li, with a 1 + 3 S[442] ground 
state, and a collection of 3 + , 2 + , and 1 + states above coming from the spin-orbit splitting 
of the 3 D[442] state. The situation is complicated by the fact that there are two linearly 
independent ways to construct an L = 2 [442] symmetry state in the p-shell. With AV18 
only, the ground state of 10 B is in fact a 1 + state, but for AV18/IL2, the spin-orbit splitting of 
the 3 D[442] states is large enough that one of the 3 + states is lowered to become the ground 
state J9j, as observed experimentally. Similar results are obtained in NCSM calculations 
using the CD-Bonn or AV8' NN potentials versus AV8' with the TM'(99) 3iV potential 
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141 ] . By comparison, 10 Be behaves more like what we expect, with a + ground 



state that is predominantly 1 S[442] symmetry in character, while the next two 2 + excited 
states are dominated by the two 1 D[442] symmetry combinations. 

The complicated spectrum for 10 B is shown in detail in Fig. El Based on GFMC calcula- 
tions, the two 3 D[442] triplets can be sorted by their quadrupole moments, Q. One triplet 
with large positive Q is widely split and contains the ground state, the second 2 + near 6 MeV 
excitation, and a predicted, but unobserved, fourth 1 + near 8 MeV. The other triplet has 
smaller negative Q and is closely spaced, starting with the second 1 + at 2 MeV excitation, 
followed by the first 2 + and second 3 + . The first 1 + is the 3 S[442] state, while the third 1 + 
around 5 MeV excitation (marked by a dash-dot line in the figure) is believed to be a 2hu 
excitation. Likewise, the second + in 10 Be near 6 MeV excitation is believed to be a 2tku 
state. These latter states will be discussed below with the unnatural-parity states. 

The A = 10 nuclei are the halfway points in the p-shell; moving further up in the shell is 
comparable to removing particles from the filled [4444] state of 16 O. The A = 11 nuclei are the 
5-hole complements of the 5- (p-shell) particle A = 9 nuclei, A = 12 nuclei are complements 
of A = 8, etc. For example, n B is the complement of 9 Be with the same allowed set of 
2S+1 L components, except that [441] symmetry becomes [443], [432] becomes [4421], [4311] 
becomes [4331], and [4221] becomes [4322]. In like manner, n Be is the complement of 9 Li, 
and n Li is the complement of 9 He. Consequently we will not give tables for these heavier 
nuclei, except for 12 C, which is of particular interest as being at the present limits of GFMC 
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, as well as being an extremely popular 



and NCSM calculations with realistic forces 
experimental target. 

Table IXlTl shows that 12 C has exactly the same 2S+l L combinations as 8 Be, with spatial 
symmetries augmented by an additional [4] in the Young diagram. The ground state will be 
a 1 5[444] + state with E OPE = 3 x (-18) = -54C three times that of 4 He. Experimentally 
12 C is 7 MeV or 8% more bound than three alphas. The distribution of the 28 pp pairs is 
such that sixteen average to give zero contribution to Eope leaving twelve ST = 10 and 
01 pairs that are equivalent to two alphas in the p-shell. In this simple model every time 
an alpha is formed in the p-shell, it effectively decouples from other nucleons in the p-shell. 
There is again a large energy gap between the first and second symmetry components in 
12 C, so we expect the ground and low-lying states to be predominantly [444] symmetry. 
By contrast, we can predict that 12 B and 12 Be should have substantial mixing of different 
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symmetries in their ground states. This knowledge is of practical benefit for the quantum 
Monte Carlo calculations, where allowing for all the possible spatial symmetries in A = 12 
nuclei is computationally prohibitive at present. 

The total energy for 30 s- and p-shell nuclei, ordered by increasing A, Z but not including 
isobaric analogs, is plotted in Fig. where experiment and our Eqpe are compared. For 
this figure, we have set the coefficient in Eq.fjHJ) to be C = 1.5 MeV. The figure shows that 
up to A ?s 9 the simple model works quite well, but then starts to underestimate the overall 
binding as A increases. Considering the necessity of including 3N forces in full GFMC 
and NCSM calculations to obtain the empirical binding energies, it is not surprising that a 
simple model based on pairwise forces will start to fail in this manner. Cohen and Kurath in 
their study of effective interactions for the p-shell U found it difficult to fit all 6 < A < 16 
nuclei at the same time, and consequently made some models to fit only A > 8 states. They 

28| that there is a gradual transition 



also found in their studies of spectroscopic factors 
from LS'-coupling to jj-coupling over the range A=9-14, and perhaps this transition is not 
unrelated to the increasing importance of 3N forces. 

As mentioned above, our simple model also predicts 10 He and 10 Li to be definitely stable, 
when they are not. In a number of other cases, the model gives identical energies for 
neighboring nuclei, such as 4 ' 5 He and 8,9 Be, and cannot predict stability one way or the 
other; this will be determined by finer details of the AW and 3N forces and kinetic energy 
considerations. Nevertheless, the simple formula reproduces the experimental trends fairly 
well. 

The model naturally indicates that total energies are close to those of summed a, t, and d 
subclusters, where applicable. In fact, the following energy relations hold for the maximally 
symmetric states with N > Z: 

E(AZ = ma) = E(AZ = ma + n) = mE a , (9) 

E(AZ = ma + 2n) = E(AZ = ma + 3n) = mE a + E 2n , (10) 

E(AZ = ma + An) = E(AZ = ma + 5n) = mE a + 2E 2n , (11) 

E(AZ = ma + 6n) = mE a + 3E 2n , (12) 

E(AZ = ma + d) = mE a + E d , (13) 

E(AZ = ma + t) = mE a + E t , (14) 

E(AZ = ma + t + n) = mE a + E t - C , (15) 
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E(AZ = ma + t + 2n) = mE a + E t + E 2n , (16) 
E(AZ = ma + t + 3n) = mE a + E t + E 2n - C , (17) 
E(AZ = ma + t + An) = mE a + E t + 2E 2n , (18) 

where m is the number of included alphas, and E 2n = —3C is the energy of a dineutron, 
which again in this simple model is equal to E^. Two non-trivial cases are Eqs. (jl5j) and 
Q17J). which in the p-shell would apply to 8,10 Li and 12 B. The model indicates there is a little 
extra binding, C more than the sum of the subclusters, on the addition of the last neutron. 
This is an accurate description of experiment for 8 Li and 12 B, but not for 10 Li. 



IV. BEYOND THE P-SHELL 



This simple model can be extended into the sd-shell, although the utility of doing so 
will continue to diminish as A increases. Counting pairs between the s-, p-, and sd-shells 
the contributions to Eope again average out so there is no net interaction between the 
shells. Then the progression from 16 to 17 0, 18 0, and 18 F nuclei is exactly analogous to 
the progression from 4 He to 5 He, 6 He and 6 Li. This is in rough accord with experiment, as 
is the prediction that the multiple-alpha nuclei will continue as Eq.Q. However, it will also 
predict that 19 F is definitely more bound than 20 O, which is not the case; among other things, 
Coulomb effects are becoming important enough that they need to be treated explicitly. 

However, the basic logic of this simple model may be applicable to sd-shell intruder 
states in the p-shell. The intruder states in A = 10 nuclei, where particle-stable intruders 
first occur, are an example. An interesting feature of the data is that the intruders in 10 Be 
are ordered starting from the most bound level as l - , 2 _ , 3 _ , and 4r . However, in 10 B the 
order is 2~, 3 _ , 4~, and then 1~. The relative ordering of the 1~ and 2~ states in these 
nuclei can be understood in the following manner. A major part of the A = 9 ground state 
is 8 Be(0 + ) plus an unpaired lp-shell nucleon (lp3/ 2 orbital in jj coupling) to which we add 
a spin-up or spin-down 2s-shell nucleon. Because these nucleons have on average no net 
interaction with the 8 Be core, their pairwise interaction should dominate. The 2~ state is a 
"stretch" state obtainable only if both spins of the pair are aligned, i.e., pure S = 1, while 
the 1~ state will have some S = pair content. In 10 Be the last pair has T = 1, so the 
2 _ state will be a 3 P pair, whereas the 1~ state will be partially a 1 S' pair, which is more 
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attractive — hence the 1~ state should be lower in the spectrum. In 10 B the last pair has 
T = 0, so the 2~ state will be a 3 S pair, whereas the 1~ state will have some admixture of 
X P, which is (much) more repulsive — hence the 1~ state will be (much) higher. Preliminary 
VMC calculations with realistic interactions successfully reproduce these level orderings and 
exhibit exactly this type of S, T pair distribution. 

One may also consider placing both last two nucleons outside the 8 Be core into the sd- 
shell, either as a dineutron pair in 10 Be or a deuteron in 10 B. In our simple model, these 
would have the same energy as the ground states, although in practice there would be some 
reduction in binding due to the greater distance from the core of these orbitals and the 
consequent overall loss of potential attraction. In actual fact, the second + in 10 Be and 
third 1 + in 10 B (shown in Fig.Elby dash-dot lines) are believed to be 2hu excitations of this 
type. These states pose an interesting challenge for both the GFMC and NCSM microscopic 
calculations. 

The present simple model provides an interesting contrast to relativistic mean-field theo- 
ries, which commonly omit the pion with the argument that its contribution will spin-isospin 
average to zero in nuclear matter; such models have been applied to nuclei as light as 16 O 
[soj ]. However, summing the expectation value of the OPE operator <7j • <Jj r» • Tj over all 
pairs we get a result that grows linearly with A for the multiple-alpha nuclei. In practice, 
the quantum Monte Carlo calculations with realistic forces find that OPE provides about 
75% of the net potential energy expectation value, although much of this comes from the 
tensor part of OPE 0. 



V. CONCLUSIONS 



We have presented a simple model for understanding the basic structure of light nuclei. 
It is based on counting the number of different S, T pairs that occur in a given nuclear 
state of specific spatial symmetry and multiplying by a numeric strength taken from one- 
pion exchange. This simple picture gives a good description for the growth of binding as A 
increases while showing saturation as the p-shell is reached. It explains the tendency of light 
nuclei to form d, t, and a subclusters, and a variety of features in the excitation spectra, 
including why for states of the same spatial symmetry, those of higher S are lower in the 
spectrum. We hope this picture provides some useful physical intuition. 
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TABLE I: Pairs and OPE 


weights 


tor A— I- 


5 nuclei. 
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5 He 
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6 
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EOPE 
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TABLE II: Ratio of kinetic to potential 


energy for A < 12 nuclei from GFMC calculations with 


the AV18/IL2 Hamiltonian. 














A Z 




Rkv 






A z 




Rkv 


2 H 




0.90 






8 Li 
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3 H 




0.84 






8 Be 
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4 He 




0.75 






9 He 
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6 He 




0.78 






9 Li 
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6 Li 




0.79 






9 Be 
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7 He 




0.79 






10 Bc 




0.77 


7 Li 




0.78 






io B 




0.77 


8 He 




0.80 






12 C 




0.77 


TABLE III: Pairs and OPE 


weights for A=6 nuclei. The Eqpe 


weight does not depend upon 


the total L value, but we enumerate all possible values for a given spin and spatial symmetry 


combination; 


the allowed T - 


- states are given under the 6 Li header, and the T = 


1 states under 


the 6 He header. 
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PP 
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01 


3 
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Pa 


6 


8 
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1 


1 


1 


EOPE 
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-3C 


+1C 
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TABLE IV: Pairs and OPE weights for A=7 nuclei. 



7 Li 7 He 
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PP 


11 




27 A 




72 


74 


72 


72 


3 


10 


3 


74 




72 


74 








01 


3 
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74 


72 






Pa 


6 


12 


Q 

o 


Q 

o 


3 


3 


3 


3 


Eqpe 
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TABLE V: Pairs and OPE weights for 8 Be. 



8 Be 
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3 


3 


72 




72 


72 
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72 




72 


72 


Pa 


6 


16 


6 


6 


6 


6 


6 
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TABLE VI: Pairs and OPE weights for 8 Li and 8 He. 









8 Li 








8 He 
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1 PDF[431] 
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72 


2 


1 
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2 


72 
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1 


72 
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72 
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1 
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Pa 


6 


6 


6 


6 


6 


6 
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Eqpe 


-10C 
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TABLE VII: Pairs and OPE weights for 9 He and 9 Li. 



9 He 9 Li 
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74 
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10 


10 


10 


10 
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TABLE VIII: Pairs and OPE weights for 9 Be. 
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10 


10 
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10 


10 
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TABLE IX: Pairs and OPE weights for 10 He and 10 Li states. 



10 He 10 Li 
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PP 
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PP 
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TABLE X: Pairs and OPE weights for 10 Be states. 
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* denotes two linearly-independent combinations 



TABLE XI: Pairs and OPE weights for 10 B states. 
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* denotes two linearly-independent combinations 



TABLE XII: Pairs and OPE weights for 12 C states. 
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Fig. 1 (Wiringa) 




FIG. 1: (Color online) Experimental spectrum for A = 6 nuclei; T = (T = 1) states are shown 
by blue (red) solid lines and breakup thresholds by black dotted lines. 
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Fig. 2 (Wiringa) 
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FIG. 2: (Color online) Experimental spectrum for A = 7 nuclei; T = 1/2 (T = 3/2) states are 
shown by blue (red) solid lines. 
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Fig. 3 (Wiringa) 




FIG. 3: (Color online) Experimental spectrum for ^4 = 8 nuclei: T = 0, 1, and 2 states are shown 
by blue, red, and green solid lines, respectively, and isospin mixed states by magenta lines. 
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Fig.4 (Wiringa) 



> 
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FIG. 4: (Color online) Experimental spectrum for ^4 = 9 nuclei: T = 1/2, 3/2, and 5/2 states 
are shown by blue, red, and green lines; solid lines denote natural-parity states and dashed lines 
unnatural-parity states. 
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Fig. 5 (Wiringa) 




FIG. 5: (Color online) Simplified experimental spectrum for A = 10 nuclei; only stable natural 
parity states are shown for 10 Be and 10 B. 



25 



Fig. 6 (Wiringa) 
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FIG. 6: (Color online) Detailed experimental spectrum for Be and B nuclei 
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Fig. 7 (Wiringa) 




FIG. 7: (Color online) Ground state binding energies for A < 16 s- and p-shell nuclei. 
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